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ABSTRACT 



Aims. Several recent studies used the hot gas fraction of galaxy clusters as a standard ruler to constrain dark energy, 
which provides competitive results compared to other techniques. This method, however, relies on the assumption that 
the baryon fraction in clusters agrees with the cosmic value f2j,/fl m , and does not differ from one system to another. 
We test this hypothesis by measuring the gas mass fraction over the entire cluster volume in a sample of local clusters. 
Methods. Combining the Sunyev-Zel'dovich thermal pressure from Planck and the X-ray gas density from ROSAT, 
we measured for the first time the average gas fraction (/gas) out to the virial radius and beyond in a large sample 
of clusters. We also obtained azimuthally-averaged measurements of the gas fraction for 18 individual systems, which 
we used to compute the scatter of / gas around the mean value at different radii and its dependence on the cluster's 
temperature. 

Results. The gas mass fraction increases with radius and reaches the cosmic baryon fraction close to ifeoo- At -R200, 
we measure / gas ,200 = 0.176 ± 0.009 (0.166 ± 0.012 for the subsample of 18 clusters in common between Planck and 
ROSAT). We find significant differences between the baryon fraction of relaxed, cool-core (CC) systems and unrelaxed, 
non-cool core (NCC) clusters in the outer regions. In average, the gas fraction in NCC clusters slightly exceeds the 
cosmic baryon fraction, while in CC systems the gas fraction converges to the expected value when accounting for the 
stellar content, without any evidence for variations from one system to another. 

Conclusions. We find that / gas estimates in NCC systems slightly disagree with the cosmic value approaching -R200. 
This result could be explained either by a violation of the assumption of hydrostatic equilibrium or by an inhomoge- 
neous distribution of the gas mass. Conversely, cool-core clusters are found to provide reliable constraints on / gas at 
overdensities A > 200, which makes them suitable for cosmological studies. 



Key words. X-rays: galaxies: clusters - 
observations 



1. Introduction 



Galaxies: clusters: general - Galaxies: clusters: intracluster medium - cosmology: 



The outskirts of galaxy clusters are the regions where 
the current activity of structure formation is taking place 
through the accretion of smaller structures onto massive 
clusters. In addition, they are the regions where the tran- 
sition between virialized gas from clusters and accret- 
ing material from the large-scale structure occurs. Non- 
ravitational effects such as turbulence (e.g., iDolag et al 



2005HVazza et all2011a) and cosmic rays (e.g. JVazza et al 
2012al iPfrommer et al.ll2007D are expected to play a more 



important role in these regions than in cluster cores, and 
the infall of smaller structures along large-scale filaments 
may cause the material in these regions to be clumpy 
dNaeai fc Laul201lT Simionescu et alJ|201lh and asymmet- 
ric (jVazza et alJ l2011bt lEckert et alJl2012j h These effects 
may bias the meas urements of cluster masses using X- 
ray and SZ proxies dfiasia et al.|[200l iNaeai et al1l2007bt 
Piffarctti & Valdarnini, l2008t ). thus setting limitations on 
the use of galaxy cluste rs as high-precision cosmological 
probes (| Allen et alJl201lh . 



Baryons in galaxy clusters, mainly in the form of hot 
X-ray emitting plasma with some contribution from stars 
in member galaxies, are subjected to gravitational and ra- 
diative processes that, in a cold dark matter (CDM) sce- 
nario, affect their distribution in the DM halo accordingly. 
By forming in the highest peaks of the primordial gravi- 
tational fluctuations, massive galaxy clusters are relatively 
well-isolated gravity-dominated structures, suggesting that 
their mass function and relative baryon budget are highly 
sensitive tests of the geometry and matter content of the 
Universe. In particular, estimates of the gas and total mass 
content of relaxed hot clusters in hydrostatic equilibrium 
can be obtained using X-ray observations, allowing one 
both to place a lower limit on the cluster baryon mass frac- 
tion, which is expected to match the cosmic value Qb/Q m 
(e.g.. I White et al. 1993]), and to constrain the cosmic dark 
energy b y imposing the ga s mass fraction as standa rd ruler 
(see e.g. lEttori et al.ll2003l 120091: lAllen et alJ[2008t ). 



In this paper, we combine the SZ pressure profiles 
from the sample of 62 galaxy clusters observed wi th the 
Planck satellite (jPlanck Collaboration et all I2012L here- 
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where /i = 0.6 is the mean molecular weight, while the gas 
mass profiles can be computed directly by integrating the 
density profiles, 
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Fig. 1. Observed gas mass fraction profiles for CC (blue) 
and NCC (red) clusters. All profiles stop at the max- 
imum detection radius (2a above the background) in 
ROSAT jPSPG. The green-shaded region shows the / gas 
profile from the average gas density and Planck pressure 
profiles. The most strongly deviating profile is that of 
A2163. 



af ter PI 2) with the R OSAT gas density profiles presented 
in lEckert et al.l (|2012t hereafter E12) to reconstruct gravi- 
tating mass profiles assuming hydrostatic equilibrium. This 
allows us to provide the first measurements of the gas frac- 
tion in cluster outskirts on relatively large samples, probing 
a volume several times larger than the typical XMM and 
Chandra limits and the entire azimuth (as opposed to most 
of the Suzaku results). We study the importance of addi- 
tional effects (non-thermal pressure support, gas clumping, 
asymmetry) on the quantities recovered by assuming that 
the gas is in hydrostatic equilibrium. Similar results are 
also presented for a subset of 18 individual clusters (6 CC 
and 12 NCC) that are in common between the samples of 
P12 and E12. In a companion paper (hereafter Paper I), we 
demonstrate the validity of the method and apply it to infer 
the average thermodynamic properties of galaxy clusters. 

Throughout the paper, we assume a ACDM cosmol- 
ogy with Sl m = 0.3, Ha = 0.7, and H a = 70 km s _1 
Mpc -1 . We also refer to the best- fit results on the ACDM 
model provided fro m WMAP-7 years (see e.g. Table 8 in 
iJarosik et al.ll201lD . which imply a cosmic baryon fraction 
^b/^m = 0.167 ± 0.011. Because our cluster sample is lo- 
cated at redshift below 0.2, the differences between these 
cosmological models affect the estimates of the gas mass 
fraction presented in this work by less than 1 per cent, 
through the variation of the angular diameter distance. 



2. Basic formalism 

For this work, the two main observables are the pressure 
from SZ measurements and the gas number density from 
X-ray observations. In this case, the hydrostatic mass is 
simply given by 



M gas (< r) = pi.m p / n gas (r')47rr' 2 dr' 



(2) 



An estimate of /gasM = M gas (< r)/Mh y d(< r) can thus be 
obtained in a straightforward way. To describe the pressure 
and the gas density, we fit the data point using general para- 
metric functions (see paper I for details). Following P12, the 
pressure is describ ed by a generalized Navarro-Frenk- Whit e 
profile (GNFW, iNavarro et aD[l997l iNaeai et al.1 [2007al) . 
The pressure gradient can be expressed analytically by dif- 
ferentiating the GNFW profile. The pressure profiles are 
taken from PI 2, for the average population as well as for 
the 18 individual clusters in common between the two sam- 
ples. 

For the density, we used a simplified version of the 
functional form introduced by IVikhlinin et al.l (|2006T ) . This 
functional form was projected along the line-of-sight and fit 
to the ROSAT emission-measure profiles. The error bud- 
get was estimated through a Monte-Carlo Markov chain 
(MCMC) method, which allows us to compute the en- 
velopes of the curves derived from the fitting procedure. 
Alternatively, we also used deprojected density profiles as 
computed in E12, and interpolated them along a common 
grid through a cubic spline method. For more details and a 
thorough validation of the method, we refer the reader to 
Paper I, also for the details of the error calculation. In Fig.[T] 
we show the gas fraction profiles derived for the 18 individ- 
ual systems and the average pressure and density profiles. 
Since we focus our work on the outer regions, the MCMC 
procedure was performed only on the data beyond 0.2i?soo, 
which causes the bumps at low radii. Beyond 0.5-Rsoo, the 
most strongly deviating profile is that of A2163, a violently 
merging cluster that is likely out of hydrostatic equilibrium 
(see Appendix C of Paper I for a discussion of this system). 

3. Results 

For a fiducial cluster with a virial temperature kT = 5 
keV, we used the universal pressure profiles from P12 to- 
gether with the median gas density profiles from El 2 to ex- 
tract gas fraction profiles using the hydrostatic equilibrium 
equation (Eq. [1]). Unlike the gas fraction profiles presented 
in Sect. 7.3 of P12, no extrapolation is needed to compute 
the gas fraction. In Fig. [2j we show the average gas frac- 
tion profile obtained by combining the average profiles from 
P12 and E12, computed through the parametric forward- 
fitting method (green-shaded area) and the geometrical 
deprojection method (red points). Excellent agreement is 
found between the results obtained between the two meth- 
ods, and hence the results presented here are stable against 
the assumed deprojection technique. The profiles obtained 
from the sample-averaged profiles are also compared to the 
median of the 18 individual objects for both deprojection 
methods. The uncertainties in the median were estimated 
using 1000 bootstrap reshuffiings of the sample (see Paper 
I). Below 1.7i?5ocQ, all methods show a remarkable level 
of agreement, demonstrating the reliability of the measure- 
ments. In the following and for the sake of simplicity we 



M hyd (< r) 



dR, 



Gnm p n gas dr 



(1) 



For a given overdensity A, 
M A /(4/37ri?i) = Ap c 



Ra is the radius for which 
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Fig. 2. (Left) Average gas fraction profile computed from the sample-average pressure and density profiles and from the 
subsample of 18 individual objects. We show the profiles obtained through parametric forward fitting with MCMC error 
envelopes for both cases (sample average, green shaded area; median of 18 objects, black circles), and by interpolating 
the deprojected density profiles (sample average, red triangles; median of 18 objects, blue circles). The error bars are 
given at the la level. (Right) Average gas fraction profiles for the CC (blue) and NCC (red) populations separately. The 
shaded areas represent the profiles obtained from the sample-averaged profiles, while the data points show the median 
of the 18 individual systems. 



focus on the results obtained with the parametric form, as 
in Paper I. However, we stress that the geometrical depro- 
jection method always leads to consistent results. 

At i?5oo, we find / gas ,500 = 0.132 ± 0.005, which 
agrees well with pre vious estimates from C handra (e.g., 
Vikh linin et al.ll2006[ ) and XMM-Newton (e.g., lEttori et all 
20101) . The measured gas frac tion reaches the cosm ic baryon 



fraction (fi 6 /fi m = 0.167, Uarosik et all 1201 if ) at R = 
1.4i?5oo and slightly exceeds it beyond this radius, although 
it remains consistent within the error bars. At R200, we 
measure 0.176 ± 0.009 (0.160 ± 0.019 in the subsample of 
18 objects) with the forward-fitting method, 0.158 ± 0.019 
(0.160 ±0.028) with the deprojection technique. Given that 
the stellar m ass contributes another 10 — 20% of the ba ryon 
fraction (e.g.. lGonzalez et al.ll2007tlLagana et al.ll2011f) . the 
gas fraction recovered assuming hydrostatic equilibrium ex- 
ceeds the expected value in cluster outskirts at the 2a 
level. Although the confidence level is low, this may in- 
dicate that either the mass recovered through hydrostatic 
equilibrium is biased low, as suggested fr om present hydro- 



dyna mical numerical si mulations (e.g., iRasia et al.l 120061 
120121 iNelson et all 120121 : iBurns et alJl2010j ). or that inho- 
moge neities in the gas distribution bias high the estimated 
M gas (jNagai fc Lavj2011t ISimionescu et al.l2011l) , or a com- 
bination of the two effects. 

Similar to the entropy (see Paper I), we observe differ- 
ences in the gas fraction profiles computed for the 6 CC and 
12 NCC clusters independently (see Fig.©. The NCC pro- 
files tend to increase more steeply in cluster outskirts com- 
pared to CC profiles. From the sample-averaged profiles, at 
i?2oo NCC clusters slightly exceed the cosmic baryon frac- 
tion (/gas, ncc — 0.169 ± 0.010 using the average profiles; 
0.178 ± 0.016 for the 12 individual NCC systems), while 
the gas fraction in CC systems is close to the expected gas 
fraction around R 200 (/ gas ,cc = 0.134±0.011; 0.143±0.015 
for the subsample of 6 CC objects). Despite its low statis- 
tical significance (less than 2a), this result likely indicates 
that hydrostatic equilibrium breaks down in the outskirts 



Table 2. Best- fitting parameters for the function given in 
Eq.Efl 



Dataset 


6gas,500 


a 


P 


X 7d.o.f 


Total 

CC 

NCC 


0.854 ±0.016 
0.757 ±0.024 
0.919 ±0.021 


— 0.20 ± 0.02 
— 0.19 ± 0.03 
-0.20 ±0.02 


0.54 ±0.05 
0.52 ±0.14 
0.46 ± 0.06 


154.4/69 
21.2/21 
100.7/45 



of NCC systems, while in CC (relaxed) clusters the behav- 
ior is closer to self-similar expectations, in agreement with 
our results for the entropy (see Paper I). 

In all cases, we can see that the enclosed gas frac- 
tion increases steadily with radius and g as temperature, as 
expec t ed from previous work (see e.g., Arnaud fc EvraTdl 



119991: lEttori fc FabiarJ 119991: IVikhlinin et al.l 120061) . To 

quantify the dependence on the gas temperature at different 
mass overdensities (200, 500, 1000, and 2500, see Fig. [3}, 
we fit a multi-linear function of the form 



fb, 



WMAP7 



5 gas,500 



( 



A 
V500 



7 keV 



(3) 



where 6 gas = / ga s//b,WMAP7 is the depletion factor. The 
best-fitting values for the parameters are provided in Table 
[2] Only the normalization &500 is found to change signifi- 
cantly between the CC and NCC populations, in agreement 
with Fig. [21 The quality of the fit, however, is found to 
change significantly when comparing the CC and NCC pop- 
ulations, with only the fit to CC objects being formally ac- 
ceptable. This indicates a significant cluster-to-cluster scat- 
ter in the latter population. 



4. Discussion 

We discuss here two important implications of our results 
on (i) the use of gas mass fraction as cosmological probe 
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Table 1. Mean / gas values and intrinsic scatter 07 (in precent) obtained from various datasets and at different overdensity 
radii. The datasets labeled "T-scaled" were rescaled by the quantity (T V i r /7 keV)^ (see Eq.[3]and Tabled]). Uncertainties 
are given at the la confidence level, while upper limits are 90%. 



Dataset 


/gas, 2500 


"7,2500 


,/gas,1000 


07,1000 


/gas, 500 


07,500 


/gas, 200 


07,200 


18 obj 


10.1 ±0.5 


1.5 ±0.4 


12.7 


± 


0.7 


2 4 +av 


14.9 


± 


1.0 


7+1. ti 

°-'-0.8 


16.6 ± 1.2 


4.2l 


1.1 
0.9 


6 CC obj 


8.7 ±0.5 


0.9 ±0.4 


10.7 


± 


0.5 


0.6 ±0.6 


12.3 


± 


0.6 


< 1.9 


14.3^ 


-1.5 
1.4 


2.8l 


1.4 
1.0 


12 NCC obj 


10.9 ±0.5 


1.2 ±0.6 


13.8 


± 


2.8 


2 - 8 ia8 


16.3 


± 


1.4 




17.8 ± 1.6 


4.7j 


1.6 
1.1 


18 obj, T-scaled 


9.9 ±0.3 


0.9 ±0.4 


12.6 


± 


0.5 


1.3 ±0.5 


14.6 


± 


0.7 


o 0+0.7 

z -°-0.6 


16.4 ± 1.1 


3.7j 


1.0 
0.8 


6 CC obj, T-scaled 


9.0 ±0.4 


< 0.8 


11.3 


± 


0.5 


< 1.0 


12.9 


± 


0.6 


< 1.3 


14.9j 


-i.fi 
1.5 


2.8l 


1.5 
1.0 


12 NCC obj, T-scaled 


10.9 ±0.5 


1 Q+0-6 

1 - v -o.5 


13.4 


± 


0.7 


1 6+ - 7 


15.6 


± 


1.0 


2.81°;? 


17.2j 


-1.5 
-1.4 


4.0l 


1.4 
1.0 




kT [keV] kT [keV] 

Fig. 3. Distribution of the gas mass fraction estimates normalized to the cosmic baryon fraction (parosik et al.|[20TTI ) as 
a function of the cluster's temperature (from E12) and core state at different overdensities. The error bars are given at 
the ltr level. The dashed lines show the best fit to the entire sample using Eq. [3] For better readability, A2163 (kT ~ 18 
keV) is not shown in the plots 



and (ii) the comparison with predictions from hydrodynam- 
ically simulated objects. 



4.1. Implications on the use of clusters for cosmology 

The use of / ya ,s as a standard ruler for cosmology 
(jEttori et all 120031: lAllen et al.l 120081 ) implies that galaxy 
clusters are representative of the cosmic baryonic budget. 
In the previous section, we quantified how the gas mass 
fraction depends on the gas temperature and matter over- 
density, providing for the first time an observationally con- 
strained calibration of the cluster baryon budget in the 
hot X-ray emitting phase. We assess here the cluster-to- 
cluster scatter of / gas at different overdensity radii. To 
compute the intrinsic scatter of a given population, we 
used a maximum-likel ihood estimator as introduced by 
iMaccacaro et al.l (|l988l) (see also Appendix A of E12). This 
method allows the mean value and the intrinsic scatter to 
be computed jointly. In Table[TJwe give the mean / gas values 
and the intrinsic scatter 07 at overdensities of 2500, 1000, 
500, and 200, for the entire population as well as for the CC 
and NCC clusters separately. We also report the scatter ob- 
tained when correcting the / gas values for their dependence 
on the cluster's virial temperature, i.e., by rescaling each 
value by the quantity (T„j r /7 keV)^ (see Eq. [3] and Table 
0). 



In Table [TJ we can see that the level of scatter increases 
with radius, from ~ 15% of the measured value at -R2500 
to ~ 25% at i?200- As expected, when a temperature scal- 
ing is performed, a lower scatter is measured in all cases. 
Interestingly, we note that the level of scatter is substan- 
tially higher for the NCC population than for the more 
relaxed, CC objects. When performing a temperature scal- 
ing, the level of intrinsic scatter is even consistent with 
in CC systems out to -R500, while at -R200 some scatter at 
the level of ~ 18% is measured, indicating that at the out- 
ermost radii even the most relaxed systems are affected by 
non- gravitational physics and/or inhomogeneous gas distri- 
butions (see E12). In the case of NCC systems, substantial 
cluster-to-cluster variations are found at all radii. 

We remark that the trend of increasing scatter with 
radius is opposite to what we would expect when prob- 
ing a larger volume, which shows that the the trend is 
likely caused by a physical effect. Both non- gravitational 
effects (turbulence, cosmic rays) and inhomogeneous gas 
densities can qualitatively explain this result. Indeed, nu- 
merical simulations show that non-thermal effe cts can be 
severa l times stronger in N CC systems (e.g., Lau et al.l 
[20091: iPfrommer etall [20071 IVazza et al.l 1201 laL I2012X 
Gas clumping and large-scale asymmetries are also ex- 
pected to be more prominent along la rge-scale filaments 
(|Nagai fc Laul 120111: IVazza et al.ll2012bl) . and thus the ef- 
fect should be stronger in merging systems. For instance, 
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if clumping biases the observed density high by ~ 20%, 
as found in Paper I for NCC clusters, the actual gas frac- 
tion for this class of objects would be consistent with the 
expected hot gas fraction. For a more detailed discussion 
of the implications of our results on the physics of cluster 
outskirts, we refer the reader to Paper I. 

To summarize, our findings (both on the global baryon 
budget and the variations in / gas from one system to an- 
other) demonstrate that the gas fraction of CC systems 
can be used efficiently as a standard ruler, while in NCC 
systems non-gravitational effects and/or inhomogeneities in 
the gas distribution introduce significant systematics in the 
computation of the gas fra ction. A very similar c onclusion 
was recently obtained by Ma hdavi et al.l ((2012) by com- 
paring X-ray and weak-lensing mass estimates. Indeed, X- 
ray mass estimates appear to be underestimated by ~ 15% 
in NCC clusters, while good agreement is found for CC 
systems. The fact that a qualitatively similar result was 
obtained using different samples and techniques thus rein- 
forces our conclusion. 

4.2. Comparison with numerical simulations 

Using the shock-capturing Eulerian adaptive mesh re- 
finem ent N-body+ga s -dyna mical adaptive refinement tree 
code, iKravtsov etaD (|2005f ) measure in the most massive 
system of M200 ~ 9.6 x 10 14 M Q a value of 6 gas that rises 
from 0.86 (at A = 2500) to 1.04 (at the virial radius) in 
the run without dissipation and from 0.56 to 0.76 when ra- 
diative cooling, star formation, metal enrichment, and stel- 
lar feedback are considered. The constraints we obtain lie 
between these predictions, suggesting that both radiative 
processes are required and the feedback provided has to be 
less incisive in its action on distributing the ICM. 

In a set of hydrodynamical simulations of massive 
galaxy clu sters performed us ing the Tr ee+SPH code 
GADGET-2, lEttori et all (|2006[ ) (see also iDe Boni et"aU 



2011) estimate at z = a & gas that increases from 0.82 
to 0.89 from A = 2500 to 200 in runs with gravitational 
heating only. On the other extreme, introducing cooling 
and star formation but no winds drastically reduces the 
expected relative gas fraction (0.29 at A = 2500; 0.59 at 
A = 200). Overall, however, none of the feedback mod- 
els explored in that work are able to reproduce the high 
6 gas value of 0.74 — 0.9 we observe at A < 500, apart 
from the cas e where the alternative arti ficial low-viscosity 
scheme from iMorris fc Monaghanl ([l997) is considered, al- 
lowing about 30 per cent of the ICM thermal energy in the 
form of turbulent motion. 

For massive halos extracted from a non-radiative 
gas-dynamical rea lization of the Millennium Simulations, 
ICrain et al.l (|2007fl measured &bar of about 0.9 (and rms of 6 
per cent) with almost no dependence on the radius between 
i?50o and -R200, implying a relative contribution from stars, 
with respect to the gas component, of 0.22 and 0.19 at the 
two radii, respectively, slightly higher than the present ob- 
ser yational constraints (e.g., < 0.17 for M 50 o > 3 x 10 14 M Q 
in Laeana et all 1201 lh . In the Millennium Gas Project, 
lYoung et al. ( 20 lit) studied the baryon distribution in 170 



objects, as SPH re-simulation of the Millennium simulation, 
with gas temperature higher than 3 keV. They considered 
models where the ICM is heated (i) solely by gravitational 
process, (ii) by pre-heating in the order of 200 keV cm 2 
occuring at z — 4 combined with cooling, (iii) by feedback 



from supernovae (SNe) and active galactic nuclei (AGN) 
as expected from the semi-analytic predictions on galaxy 
formation but not including radiative cooling. For the most 
massive systems, they measured 6 gas in the range 0.88-0.91, 
0.57-0.83, 0.70-0.82 between A = 2500 and 200 for cases 
(i), (ii), (iii), respectively, suggesting that (a) some cooling- 
is needed to lower the gas fraction predicted in the inner 
part and (b) the implemented action of either pre-heating 
or feedback from SNe and AGN is probably too strong to 
fully explain the observed gas fraction in the clusters' out- 
skirts. 



5. Conclusion 

Combining X-ray (ROSAT/PSPC) and SZ (Planck) data, 
we analyzed the properties of the ICM at large radii using 
both sample- averaged profiles and the data for 18 individ- 
ual systems. This allowed us to measure for the first time 
hydrostatic mass and gas fraction profiles out to the virial 
radius in a relatively large sample of objects. The results 
on the thermodynamic properties of galaxy clusters close to 
the virial radius are presented in a companion paper. Our 
results can be summarized as follows: 

— The gas fraction reconstructed through the hydrostatic 
equilibrium equation increases with radius and reaches 
the value of / gas ,2oo = 0.176 ± 0.009 around i?2oo- 
Different input data and deprojection techniques yield 
results that excellently agree. This value is consis- 
tent with the cosm ic baryon fraction f2{,/f2 m = 0.167 
(|Jarosik et al.l [20TTt ) . Given that the stellar content is 
expected to account for 10 — 20% of the baryons in 
galaxy clusters, the measured gas fraction slightly ex- 
ceeds the expected value. A 15% bias either in raising 
the hydrostatic masses, as expected from numerical sim- 
ulations, or in lowering the gas mass allows us to recon- 
cile our measurement with the expected gas fraction. 

— We observe differences between the gas fraction mea- 
sured in relaxed, CC systems compared to dynamically 
active, NCC systems. In NCC clusters, the gas frac- 
tion exceeds the cosmic baryon fraction significantly, 
indicating that hydrostatic equilibrium breaks down in 
the outer regions of perturbed systems, and/or that the 
gas is distributed in an inhomogeneous way. Conversely, 
in CC systems the gas fraction converges to the ex- 
pected hot gas fraction, confirming that they correspond 
to more relaxed systems, which can be used more effi- 
ciently for cosmological purposes. 

— From the 18 objects comprising our sample, we mea- 
sured the scatter around the average / gas value using a 
maximum-likelihood estimator. We find that the scat- 
ter increases with radius, from 15% at -R2500 to 25% at 
-^200- In contrast to NCC systems, in which we measure 
substantial cluster-to-cluster variations at all radii, CC 
systems exhibit a negligible level of scatter out to -R500, 
confirming that they are better suited for cosmological 
studies. 

— We quantified the / gas dependence upon gas tem- 
perature and matter overdensity with the relation 
/ g as//b,WMAP7 = ^500 (A/500)" (T gas /7 keV) 13 and 
measured 6500 = 0.76 ± 0.02 in CC systems and 0.92 ± 
0.02 in NCC clusters, a = -0.2 and /3 = 0.5. The gas 
bias (or depletion) factor 6500, and its variation with 
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temperature and radii, disagrees with the results ob- 
tained for clusters extracted from recent hydrodynami- 
cal simulations, implying that the action of cooling and 
AGN and/or SN feedback is indeed needed but has to 
be regulated differently from what has been used so far. 

References 

Allen, S. W., Evrard, A. E., & Mantz, A. B. 2011, ARA&A, 49, 409 
Allen, S. W., Rapetti, D. A., Schmidt, R. W., et al. 2008, MNRAS, 
383, 879 

Arnaud, M. & Evrard, A. E. 1999, MNRAS, 305, 631 
Burns, J. O., Skillman, S. W., & O'Shea, B. W. 2010, ApJ, 721, 1105 
Crain, R. A., Eke, V. R., Frenk, C. S., et al. 2007, MNRAS, 377, 41 
De Boni, C, Dolag, K., Ettori, S., et al. 2011, MNRAS, 415, 2758 
Dolag, K., Vazza, F., Brunetti, G., & Tormen, G. 2005, MNRAS, 364, 
753 

Eckert, D., Vazza, F., Ettori, S., et al. 2012, A&A, 541, A57 
Ettori, S., Dolag, K., Borgani, S., & Murante, G. 2006, MNRAS, 365, 
1021 

Ettori, S. & Fabian, A. C. 1999, MNRAS, 305, 834 
Ettori, S., Gastaldello, F., Leccardi, A., et al. 2010, A&A, 524, A68 
Ettori, S., Morandi, A., Tozzi, P., et al. 2009, A&A, 501, 61 
Ettori, S., Tozzi, P., & Rosati, P. 2003, A&A, 398, 879 
Gonzalez, A. H., Zaritsky, D., & Zabludoff, A. I. 2007, ApJ, 666, 147 
Jarosik, N., Bennett, C. L., Dunkley, J., et al. 2011, ApJS, 192, 14 
Kravtsov, A. V., Nagai, D., & Vikhlinin, A. A. 2005, ApJ, 625, 588 
Lagana, T. F., Zhang, Y.-Y., Reiprich, T. H., & Schneider, P. 2011, 
ApJ, 743, 13 

Lau, E. T., Kravtsov, A. V., & Nagai, D. 2009, ApJ, 705, 1129 
Maccacaro, T., Gioia, I. M., Wolter, A., Zamorani, G., & Stocke, J. T. 

1988, ApJ, 326, 680 
Mahdavi, A., Hoekstra, H., Babul, A., et al. 2012, ArXiv e-prints 
Morris, J. P. & Monaghan, J. J. 1997, Journal of Computational 

Physics, 136, 41 

Nagai, D., Kravtsov, A. V., & Vikhlinin, A. 2007a, ApJ, 668, 1 
Nagai, D. & Lau, E. T. 2011, ApJ, 731, L10+ 
Nagai, D., Vikhlinin, A., & Kravtsov, A. V. 2007b, ApJ, 655, 98 
Navarro, J. F., Frenk, C. S., & White, S. D. M. 1997, ApJ, 490, 493 
Nelson, K., Rudd, D. H., Shaw, L., & Nagai, D. 2012, ApJ, 751, 121 
Pfrommer, C, EnBlin, T. A., Springel, V., Jubelgas, M., & Dolag, K. 

2007, MNRAS, 378, 385 
Piffaretti, R. & Valdarnini, R. 2008, A&A, 491, 71 
Planck Collaboration, Ade, P. A. R., Aghanim, N, et al. 2012, ArXiv 

e-prints 

Rasia, E., Ettori, S., Moscardini, L., ct al. 2006, MNRAS, 369, 2013 
Rasia, E., Meneghetti, M., Martino, R., ct al. 2012, New Journal of 

Physics, 14, 055018 
Rasia, E., Tormen, G., & Moscardini, L. 2004, MNRAS, 351, 237 
Simionescu, A., Allen, S. W., Mantz, A., et al. 2011, Science, 331, 

1576 

Vazza, F., Briiggcn, M., Ghcller, C, & Brunetti, G. 2012a, MNRAS, 
421, 3375 

Vazza, F., Brunetti, G., Gheller, C, Brunino, R., & Briiggen, M. 

2011a, A&A, 529, A17+ 
Vazza, F., Eckert, D., Simionescu, A., Brueggen, M., & Ettori, S. 

2012b, ArXiv e-prints 
Vazza, F., Roncarelli, M., Ettori, S., & Dolag, K. 2011b, MNRAS, 

413, 2305 

Vikhlinin, A., Kravtsov, A., Forman, W., et al. 2006, ApJ, 640, 691 
White, S. D. M., Navarro, J. F., Evrard, A. E., & Frenk, C. S. 1993, 
Nature, 366, 429 

Young, O. E., Thomas, P. A., Short, C. J., & Pearce, F. 2011, MNRAS, 
413, 691 



6 



